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ABSTRACT 

We present the analysis of the X-ray Multi-Mirror Mission {XMM-Newton) European Photon 
Imaging Camera (EPIC) data of the Galactic supernova remnant (SNR) CTB 109 (G109. 1-1.0). 
CTB 109 is associated with the anomalous X-ray pulsar (AXP) IE 2259-1-586 and has an unusual 
semi-circular morphology in both the X-ray and the radio, and an extended X-ray bright interior 
region known as the 'Lobe'. The deep EPIC mosaic image of the remnant shows no emission 
towards the west where a giant molecular cloud complex is located. No morphological connection 
between the Lobe and the AXP is found. We find remarkably little spectral variation across the 
remnant given the large intensity variations. All spectra of the shell and the Lobe are well fitted 
by a single-temperature non-equilibrium ionization model for a coUisional plasma with solar 
abundances (fcT k, 0.5 - 0.7 keV, r = / n^dt « 1 - 4 x 10" s cm'^, « 5 - 7 x lO^i cm^^)^ 
There is no indication of nonthermal emission in the Lobe or the shell. We conclude that the 
Lobe originated from an interaction of the SNR shock wave with an interstellar cloud. Applying 
the Sedov solution for the undisturbed eastern part of the SNR, and assuming full equilibration 
between the electrons and ions behind the shock front, the SNR shock velocity is derived as 
Vs = 720 ± 60 km s~^, the remnant age as t — (8.8 ± 0.9) x 10'^ d^ yr, the initial energy as 
Eq = (7.4 ± 2.9) X 10^" d\-^ ergs, and the pre-shock density of the nuclei in the ambient medium 
as no = (0.16 ±0.02) dg^'^ cm~'^, at an assumed distance oi D = 3.0 d^ kpc. Assuming CTB 109 
and IE 2259+586 are associated, these values constrain the age and the environment of the 
progenitor of the SNR and the pulsar. 

Subject headings: shock waves — supernova remnants — X-rays: individual (CTB 109) — ISM: individual 
(CTB 109) 
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1. Introduction 

The Galactic supernova remnant (SNR) CTB 109 
(G109. 2-1.0) was discovered in X-rays with Ein- 
stein by Gregory & Fahlman (1980). In the radio 
band, it was identified as an SNR by Hughes et al. 
(1981) in the Galactic plane survey at A49 cm 
with the Westcrbork Synthesis Radio Telescope. 
CTB 109 is the host remnant of the anomalous 
X-ray pulsar (AXP) IE 2259+586 with a spin pe- 
riod of P = 6.98 s (Fahlman & Gregory 1981, 
1983) and a spin-down rate of 3 — 6 x 10"^'^ s s~^ 
(Iwasawa et al. 1992). In 2002, an outburst of 
IE 2259-1-586 occurred with bursts similar to Soft 
Gamma Repeaters and with a sudden spin up of 
6u/u = 4 X 10"^ (Kaspi et al. 2003; Gavriil et al. 
2003; Woods et al. 2004). 

The distance to the SNR has been estimated 
from various observations as 3 - 6 kpc (Sofue et al. 
1983; Hughes et al. 1984; Kothes et al. 2002, and 
references therein). A distance of 3.0 ±0.5 kpc has 
been derived by Kothes et al. (2002) by measuring 
the spectroscopic distances and radial velocities of 
H II regions and comparing those values to the ra- 
dial velocity of CTB 109 measured by Tatematsu 
et al. (1990). The remnant is embedded in a large 
complex of H II regions which extends over 400 pc 
along the Galactic plane. 

At radio frequencies, Hughes et al. (1984) have 
derived a spectral index of a = 0.50 ± 0.04, for 
flux density Si, oc and frequency v, which 

does not vary across the remnant. The remnant 
shell is incomplete in the west, both in radio and 
X-rays (see Fig. 1). Bright spots arc found aromid 
the rim, but the features do not correlate in the 
two bands. No radio point source is found at the 
position of the X-ray pulsar, which is displaced by 
3.'6 from the geometrical center of the radio shell. 

CTB 109 is located near a giant molecular cloud 
(GMC) complex (Israel 1980) which contains five 
Hll regions. Heydari-Malayeri et al. (1981) and 
Tatematsu et al. (1985) have suggested that this 
GMC complex is associated with the SNR based 
on CO and ^^CO observations. Data from the Po- 
sition Sensitive Proportional Counter (PSPC) of 
the Rontgensatellit (ROSAT) confirm that there 
is no X-ray emission in the west (Rho & Petre 
1997, hereafter RP97). The SNR shell does not 
extend further to the west behind the GMC in 
either radio or X-rays; because absorption is neg- 



ligible in the radio continuum, this implies that 
the semi-circular shape of the remnant is not due 
to absorption. More likely, the SNR shock has 
been stopped by the GMC complex on one side. 
To the east, the interstellar medium (ISM) density 
appears to be lower, as indicated in the H i map of 
Kothes et al. (2002). In CO (J = 1 - 0) and ^^CO 
( J = 1 — 0) observations, Tatematsu et al. (1987) 
have found an arm-like CO ridge ('CO arm') which 
is anti-correlated with the X-ray emission. 

CTB 109 has a strikingly bright region in X- 
rays in the interior, the 'X-ray blob' or the 'Lobe'. 
Gregory & Fahlman (1983) have suggested that 
the AXP emits a jet and that the bright Lobe is 
material excited by a jet or the jet itself. However, 
Rho & Petre (1993, 1997) have shown that the 
ROSAT PSPC spectrum of the Lobe is thermal in 
origin and shows no evidence of synchrotron emis- 
sion, arguing against a jet interpretation. Hur- 
ford & Fesen (1995) have found no morpholog- 
ical evidence to support a pulsar jet origin for 
the Lobe in the ROSAT High Resolution Imager 
(HRI) data and conclude that the bright X-ray 
emission from the Lobe is most likely caused by 
a density enhancement in material stripped from 
the cloud complex. Finally, the Chandra image 
(Patel et al. 2001) reveals no morphological con- 
nection between the Lobe and the AXP. Rho & 
Petre (1997) show that for parts of the shell, the X- 
ray spectrum is well fitted by a single-component 
thermal model, whereas for the Lobe and the 
northern and the southern part of the shell, two 
thermal components are required. A BeppoSAX 
Low-Energy Concentrator Spectrometer (LEGS) 
spectrum of the whole SNR excluding the pul- 
sar and the Lobe has been analyzed by Parmar 
et al. (1998). The parameters of the best fit to 
the LEGS data with a non-equilibrium model are 
kT = 0.95tg:^^ keV and the ionization timescale 
net = 3.8tli X 10^^ cm-3 s. 

Surveys of Galactic SNRs have been performed 
with the Infrared Astronomical Satellite (IRAS), 
revealing infrared emission from CTB 109 (Arendt 
1989; Sakcn et al. 1992). The data in the wave- 
length bands of 12, 25, 60, and 100 i^m show a 
dense dust cloud which is correlated with bright 
H II regions, located in the south and west of the 
remnant (Goe et al. 1989). Furthermore, an en- 
hancement in infrared emission is found between 
the Lobe and the north-east shell. 
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Optical emission is detected mainly near the 
bright northeastern and southern radio emission 
(Hughes et al. 1981). Spectra of filamentary struc- 
tures in the south confirm their shock origin (Blair 
& Kirshner 1981). Ha and [Oili] images of Fesen 
& Hurford (1995) reveal previously unknown fila- 
ments, e.g. faint filaments near the projected cen- 
ter of the remnant, along the southwestern edge 
of the bright Lobe. Prom the filament spectra, 
they derive a shock velocity of ^ 100 km s"""^ and 
pre-shock cloud densities of 5 - 20 cm~^. 

We have undertaken observations of CTB 109 
with the XMM-Newton observatory to provide the 
deepest X-ray image of this SNR and to provide 
the highest quality X-ray spectral data to date. 
The primary objectives of our analysis arc to de- 
termine the evolutionary parameters of the rem- 
nant assuming it is in the Sedov phase and to de- 
termine the nature of the Lobe emission. The deep 
X-ray images are used to study the morphology of 
this remnant in comparison to the data at other 
wavelengths and to search for any morphological 
connection between the pulsar and the Lobe. The 
high quality spectral data are used to determine 
the evolutionary parameters of the remnant (age, 
shock velocity, ambient density) and to search for 
spectral variations within the remnant. The evo- 
lutionary parameters determined from the X-ray 
spectral fits provide useful constraints not only on 
the age and the environment in which the remnant 
formed but also for the pulsar. This is particu- 
larly useful for the AXP since there is no other 
reliable technique for estimating the age of this 
object (Kaspi et al. 2003; Woods ct al. 2004). The 
detailed spectral fits of the shell and the Lobe pro- 
vide useful insights into the nature of the Lobe 
emission. In particular, the high-quality spectra 
of the Lobe are used to determine if the abun- 
dances are enhanced (evidence of eject a), if there 
is a nonthermal component, and if the tempera- 
ture and ionization state of the plasma is similar 
to that of the shell. 

2. XMM-Newton Data 

We have observed CTB 109 with three point- 
ings of the X-ray Multi-Mirror Mission (XMM- 
Newton, Janscn et al. 2001; Aschenbach et al. 
2000) in the AOl period covering the southern 
(hereafter: pointing S), northern (pointing N), and 



eastern parts (pointing E). We analyze data from 
the European Photon Imaging Cameras (EPIC) 
PN (Striider et al. 2001), EPIC MOSl, and EPIC 
M0S2 (Turner et al. 2001). The exposure times, 
after removing flares and periods of high back- 
ground, are ~10 ks. Furthermore, we have access 
to observation 0038140101 (PI: V. Kaspi, here- 
after pointing PI), with a net integration time 
of 34 ks for MOSl/2. We also analyze the tar- 
get of opportunity (TOO) dataset of the pulsar 
IE 2259-F586 in the XMM-Newton Archive (Obs. 
ID 0155350301, hereafter P2) with 17 ks of M0S2 
data in full frame mode. The MOSl was in partial 
window mode, covering only the pulsar. In obser- 
vations PI and P2, the PN was in small window 
mode. Therefore, wc do not consider the PN data 
of the observation PI and the PN and MOSl data 
of observation P2 in our analysis. Observation de- 
tails of the datascts used arc listed in Table 1. 
Starting from the observational data files (ODFs), 
the data are processed with XMM-Newton Science 
Analysis System (XMMSAS) version 5.4.1. For 
EPIC PN, only single and double pattern events 
are used, whereas for the MOSl and 2, singles to 
quadruples are selected. 

2.1. Intensity Map 

The SNR has been fully covered in five point- 
ings with different directions (three for the shell 
and two centered on the pulsar). We therefore 
merge the images from the different pointings to 
obtain complete exposure corrected images of the 
remnant using the procedure emosaic. For all the 
pointings, events are selected in the energy bands: 
0.3 - 0.9 keV, 0.9 1.5 keV, 1.5 - 4.0 keV, and 
4.0 keV - 12.0 keV. Images are created with a bin 
size of 4" and smoothed with circular Gaussian 
function that has a sigma of 3 pixels in both di- 
mensions RA and Dec. In the 4.0 keV - 12.0 keV 
images, only significant emission from the pulsar 
is seen, but not from the rest of the SNR. 

Figure 1 shows a mosaic image in the energy 
band of 0.3 - 4.0 keV in false color. The rem- 
nant's X-ray emission fades off to the west, lacking 
a clearly defined edge. In the interior, the Lobe is 
brighter than any part of the shell and there are 
dark regions north and south of the Lobe. The 
Lobe emission is well separated from the pulsar. 
We estimate the distance of the shell from the pul- 
sar in the mosaic image by fitting a circle centered 
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on the pulsar to each sector and obtain: 18'.5 for 
the eastern shell, 17.'3 for the southern shell, and 
16^6 for the northern shell. Thus, the eastern shell 
has the largest projected distance from the pulsar. 

2.2. RGB Composite Image 

As no emission from the remnant is found above 

4 keV except from the pulsar, we create images in 
the energy bands below 4 keV in order to identify 
possible spectral variations within the remnant. 
The images in the bands red (R) = 0.3 - 0.9 keV, 
green (G) = 0.9 - 1.5 keV, and blue (B) = 1.5 
- 4.0 keV are combined into an RGB mosaic im- 
age (Fig. 2). Wc then use this image to guide our 
selection of regions to be examined more closely 
through detailed spectral fitting. The selected re- 
gions are shown in Figure 3. 

The shell is slightly redder than the interior and 
has some red clumps. There is a small sector of the 
shell with lower emission in the northeast, coin- 
ciding with a bright region in the radio continuum 
(Kothes et al. 2002). Next to this region, there are 
red knots in the RGB composite image of the X- 
ray shell (regions 5 and 6 in Fig. 3). Furthermore, 
there is a larger red spot in the south (region 15). 
There arc also color variations within the Lobe, 
the outskirts being redder than the central part of 
the Lobe. 

In addition to the very bright pulsar, fainter 
point sources are found in and around the rem- 
nant. Two point sources located southeast of the 
remnant appear red in the RGB composite image. 
The one further east at RA = 23''03'"19.6^ Dec = 
-F58°45'29" (J2000.0) is a star of the spectral type 
G5.5:V (USNO-A2.0 1425-14515707, Monet et al. 
1998) and is coincident with source No 3 in RP97. 
The other source at RA = 23''02'"43.7^ Dec = 
-^58°37'46" (J2000.0) is coincident with source No 
4 in RP97; it is also coincident with a point source 
in the DSS image. Presumably, it is also a Galac- 
tic star. The point sources close to the pulsar 
(in projection) have blue or violet colors. The 
source at RA = 23'»00'"33.3^ Dec = -F58°52'45" 
(J2000.0) located west of the pulsar is USNO- 
A2.0 1425-14436845 (No 1 in RP97), a star of 
the spectral type G9.5:III-IV. The source at RA = 
23''00'"43.3", Dec = +58°50'28" (J2000.0) south 
of USNO-A2.0 1425-14436845, is only known as an 
X-ray source (No 2 in RP97). There is an optical 
point source on the DSS image at the X-ray posi- 



tion. Therefore, these hard point sources seem to 
be Galactic stars as well. Enhanced column den- 
sity in the foreground or high intrinsic absorption 
seems to be responsible for the hard spectrum. We 
extracted spectra for these point sources, but the 
low number of counts does not allow detailed spec- 
tral fitting. 

2.3. Background 

Since CTB 109 is a relatively faint diffuse ob- 
ject, a good knowledge of the background is cru- 
cial for the spectral analysis. We studied the 
background properties of the EPIC cameras based 
on the treatments of Lumb (2001), Lumb et al. 
(2002), and Read & Ponman (2003). 

In the case of a point source a local background 
extracted from another region of the same data 
as the source can be used. For extended objects, 
it is not possible to obtain background data at 
small distances from the source, i.e. a few arcsec- 
onds away from the source like for point or point- 
like sources. Instead one has to extract the back- 
groimd data far from the target, from a part of the 
detector that is not illuminated by source emis- 
sion. In this case the distance between the source 
and the background region is large (few arcmin- 
utes) and the difference of the detector background 
at the two positions could become significant (see 
next paragraph). In the case of a diffuse source 
that fills the entire field of view, it is even impos- 
sible to estimate the local background from the 
same data. If one selects a background region far 
from the source itself, many effects can produce an 
inappropriate background due to the different chip 
position: First, the effective area of the mirrors 
depends on the off-axis angle. Photons are sub- 
ject to vignetting, but particles are not. Second, 
high-energy particles which interact with material 
surrounding the detector produce fluorescence, the 
magnitude of this fluorescence emission can vary 
significantly with position on the detector, espe- 
cially for the PN detector. Third, the spectral 
response depends on the position on the detector, 
especially for the PN. 

The internal 'quiescent' background of the 
EPIC cameras consists of fluorescence lines and 
events produc;cd by charged particles. The flu- 
orescence depends on the camera body materials 
and the incident flux of high energy charged parti- 
cles. Therefore, the fluorescent component varies 
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with position on the detector and with time. In 
addition, the high energy charged particles can 
produce background events directly in the detec- 
tor. Finally, low energy protons can reflect off 
of the mirror surfaces and reach the detectors in 
the focal plane. The background produced by low 
energy protons is highly variable in time for the 
XMM-Newton instruments. 

There are two types of background template 
files available that can be used for the analysis 
of EPIC data: 'Closed' data and blank sky data. 
'Closed' data have been obtained with the filter 
wheel in the 'closed" position; these data include 
only the internal detector background. However, 
this instrument configuration is not representative 
of the usual observation where the CCDs are ex- 
posed to photons and low energy charged parti- 
cles. For the blank sky data, datasets for different 
pointings have been merged for each EPIC cam- 
era after eliminating sources from the data. These 
datasets comprise the detector background and an 
average cosmic X-ray background. Since the back- 
ground component due to high energy charged 
particles is included in these data at an unknown 
level, we need to model it using the spectral shape 
determined from the closed data. 

2.4. Extraction of EPIC Spectra 

2.4- 1- Circular Regions 

Although the cosmic X-ray background is a sig- 
nificant uncertainty in the blank sky data, we de- 
cide to use the blank sky instead of the 'closed' 
data, because the latter have poor photon statis- 
tics and are therefore not suitable for extracting 
small regions. We use the blank sky data of the 
background analysis group in Birmingham, UK 
because they have data for both the medium and 
the thin filters, as well as for the full frame and 
extended full frame PN modes. The Birmingham 
group have also developed scripts that can be used 
for the background analysis. The script skycast 
converts the detector coordinates of the template 
background datasets into sky coordinates using 
the pointing direction of the observation to be 
analyzed and adds sky coordinates to the back- 
ground data. The script createspectra extracts 
spectra in specified circular (or annular) regions 

See http://www.sr. bharn. ac.uk/xrnm3/. 



from the data and background-file. Ancillary re- 
sponse files (ARFs) and redistribution matrix files 
(RMFs) are produced for the corresponding detec- 
tor regions using the XMMSAS commands arf gen 
and rmf gen. Since the source spectrum and asso- 
ciated blank sky background spectrum are created 
at the same position of the detector, the ARF and 
RMF are correct for both source and background. 

Spectra are analyzed with XSPEC Ver. 11.3. 
We read in the source and the blank sky back- 
ground spectra as two datasets: source (1), back- 
ground (2). We model the background spectrum 
as a combination of a powerlaw component and 
two Raymond & Smith (Raymond & Smith 1977, 
hereafter RS) components with temperatures kT 
= 0.08 keV and 0.2 keV (Lumb 2001; Lumb et al. 
2002). Zero- width Gaussians are used to repro- 
duce the fluorescence lines of the detector. The 
detector background is modeled with a powerlaw 
component and Gaussians using the XSPEC op- 
tion '/b' to suppress the application of the ARF. 
An additional thermal component is used to model 
the SNR emission (set to zero for the background 
dataset). In Figure 4, the EPIC PN and MOSl/2 
spectra of region 1 (see Fig. 3) extracted from the 
E pointing are displayed together with the corre- 
sponding blank sky background spectra. We show 
the spectrum from region 1 since this is the faintest 
part of the Lobe, and hence the most sensitive to 
the background subtraction. 

In order to convince ourselves that the back- 
ground in the source spectrum can be well es- 
timated from the blank sky data, we compare 
the spectra extracted as source (from observa- 
tions) and background (from blank sky data) in 
a region outside the SNR. Both for PN and for 
MO S 1/2, the off-remnant spectrum from observa- 
tions matches the blank sky spectrum very well 
in the softer (< 0.6 keV) and harder (> 1.5 keV) 
bands as well as for the ffuorescence lines. There 
is a small excess in the off-remnant spectrum be- 
tween ^ 0.6 and 1.5 kcV. This faint emission is 
well fitted with a thermal model component and 
is likely Galactic diffuse X-ray emission. Since 
the remnant spectrum dominates the background 
spectrum in the 0.6 - 1.5 keV band (see Fig. 4), 
using the blank sky data as background intro- 
duces only a small systematic uncertainty in the 
derived parameters for the remnant spectra. We 
verify this by modeling the spectra of regions in 
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the remnant shell including an additional thermal 
component for the local diffuse X-ray background, 
as derived from the off-remnant spectrum. The 
parameters of the model component for the rem- 
nant emission (foreground absorption, tempera- 
ture, and ionization timescale, see § 3) are consis- 
tent with results from the fit using only the blank 
sky data as background. The abundances of the 
elements with emission lines in this energy band 
(O, Ne, and Mg) also yield fit results that agree 
well with the results without the additional com- 
ponent, i.e. solar abundances. Consequently, the 
use of the blank sky data as background should 
have little effect on our determination of the neu- 
tral hydrogen column density and on the temper- 
ature of the thermal component. 

2.4-2. Arbitrary Extraction Region Shape 

Since the existing version of arf gen only com- 
putes ARF files for circular and annular extrac- 
tion regions, our analysis has been largely limited 
to circular regions. However, the XMMSAS com- 
mand evigweight calculates the vignetting cor- 
rection for each event and thus one can use the on- 
axis ARF. In order to subtract the non- vignetted 
events from the observation data, it is better to 
use the 'closed' data for the background, as this 
dataset only consists of the non- vignetted compo- 
nent, which is the internal detector background, 
and thus includes no additional external back- 
ground emission. Compared to the blank sky data, 
the 'closed' data have far poorer statistics. There- 
fore, we use this method only for non-circular ex- 
traction regions. 

After computing the weight factor for the vi- 
gnetting correction for each event, spectra are ex- 
tracted from the new event files and RMFs are 
created with rmfgen. The source data still in- 
clude the non-vignetted events which have not 
passed through the X-ray telescope, but now have 
a weighting factor applied. This misapplication 
of the weighting factor to non-vignetted events 
can be compensated for by applying the same 
weights to the 'closed' background data. Coor- 
dinates of the 'closed' data are converted into the 
corresponding sky position using skycast. The vi- 
gnetting correction weights arc computed for the 
'closed' data. Finally, spectra are extracted in 
the same regions as for the observed source data. 
Again the background spectrum is used as second 



dataset and modeled as in § 2.4.1. 

2.4.3. Comparison of the Two Methods 

In order to check for consistency between the 
two methods for creating spectral data, spectra of 
the Lobe regions (see § 3.4), extracted from the 
pointing E data, are analyzed using both meth- 
ods and compared. For the first method, we 
use the blank sky data as background data and 
create the ARF with arfgen. For the second 
method, the event files are vignetting-corrcctcd 
with evigweight and on-axis ARF and 'closed' 
background data are used. All the fitted param- 
eters derived from the two methods agree within 
90% confidence limits. Therefore, we are confi- 
dent that either method may be used without in- 
troducing systematic errors which dominate over 
the statistical errors. 

3. The Spectral Analysis 

The XMM-Newton EPIC data allow us to per- 
form spatially resolved spectral analysis of the 
SNR. Extraction regions are shown in Figure 3 
on the 0.3 4.0 kcV mosaic image. Datasets and 
extraction regions are selected as follows: 

• As we have shown in § 2.1, the eastern part 
of the shell is farthest from the pulsar. Since 
it is located opposite to the CMC complex 
in the west, the eastern shell segment most 
likely indicates a blast wave propagating into 
a lower-density ISM (Kotlics et al. 2002). 
For the spectral analysis, we cover the east- 
ern shell with three circular regions (7 - 9) 
and also arc regions for comparison. In addi- 
tion, we extract spectra for regions next to 
the radio bright knot in the northeast (re- 
gions 5 and 6). Dataset E is used for all 
these analyses. 

• The shell is brightest in the northern part; it 
has been covered completely by pointing N. 
We select circular regions from this dataset 
centered on the brighter knots (regions 12 
and 13) and also analyze interior regions 
close to the northern shell for comparison 
(regions 11 and 14). 

• In the southern part of the shell which does 
not have as clearly a defined edge as the east- 
ern or the northern portions of the shell, we 
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study the EPIC spectra (using observation 
S) for two regions (15, 16) with soft emis- 
sion (see RGB composite image, Fig. 2). 

• In CO observations, an enhancement of CO 

gas has been found which extends from the 
GMC in the west across the northern part of 
the remnant (Tatematsu et al. 1987). This 
dense material is called the CO arm. Its po- 
sition and shape coincides well with an X-ray 
faint interior region, north of the pulsar and 
the Lobe (see Fig. 1). Spectra are analyzed 
for a large circular region (10) as well as for a 
non-circular region which follows the shape 
of the very low emission region. Data from 
observations N and PI are used. 

• The interior of the remnant has remarkably 
low surface brightness in the south ('Dark 
Interior'). This region is analyzed using the 
PI data which have the longest exposure 
time. 

• We extract spectra for four regions of the 
Lobe (regions 1 - 4 in Fig. 3) from the ob- 
servations E, PI, and P2, in order to look 

for spectral variations. 

All the EPIC spectra are analyzed in XSPEC us- 
ing an energy range of 0.3 - 10.0 keV. The spectra 
arc grouped with a minimum of 20 counts per bin 
and the statistic is used. The errors quoted are 
90% confidence limits. As can be seen in Figure 
4, the emission above ^ 4 keV is mainly back- 
ground emission and is crucial for estimating the 
background. Although we select the energy range 
of 0.3 - 10.0 keV for spectral fitting, we will only 
show the spectrum below 4 keV in the following 
figures in order to clearly display the SNR emis- 
sion. 

For analyzing the XMM-Newton EPIC spectra 
we use the following models: VMEKAL, which de- 
scribes an emission spectrum for hot ionized gas 
in collisional ionization equilibrium with variable 
abundances (Meweet al. 1985, 1986; Kaastra 1992; 
Liedahl et al. 1995), and VNEI, which is a model 
for a non-equilibrium ionization (NEI) collisional 
plasma, again with variable abundances (Hamil- 
ton et al. 1983; Borkowski et al. 1994; Liedahl et al. 
1995; Borkowski et al. 2001). The VNEI model 
includes a parameter for the ionization timescale, 
T = J riedt, an integral over the time since the gas 



was shocked. In addition, wc modeled some re- 
gions in the shell with the model VPSHOCK for a 
constant temperature, plane-parallel shock which 
includes a range of ionization timcscalcs, parame- 
terized by the lower limit and the upper limit of 
the T range (Borkowski et al. 2001, and references 
therein) . 

The XMM-Newton EPIC spectra of region 3 
in the Lobe extracted from the pointing E data 
(with VNEI model fits) are presented in Figure 
5. Emission line features corresponding to the 
Mgxi triplet and the Sixiii triplet are easily vis- 
ible. Moreover, the Nex Lya line can be clearly 
seen in the MOS data, as the energy resolution of 
MOS is slightly better than that of PN. 

We first fit the spectra with free abundances. 
For all the spectra, the abundances for all the el- 
ements are consistent with solar abundances. We 
thereupon fix the abundances to solar values in or- 
der to reduce the number of free parameters and 
hopefully to better constrain the remaining free 
parameters. Only Mg and Si are still kept free, 
because these abundances have appeared to differ 
slightly from solar values for some regions and the 
line features are evident in the EPIC spectra. We 
allowed the Si and Mg abundances to vary inde- 
pendently. For other elements like Ne or S, it is too 
difficult to distinguish between line and continuum 
emission in the corresponding energy intervals. 

All the EPIC spectra are fitted better by a 
single-temperature VNEI model than by a colli- 
sional ionization equilibrium VMEKAL model. In 
the following, we present only the VNEI results 
for all the regions, and discuss the VPSHOCK re- 
sults for shell spectra where we expect the least 
disturbed plane-parallel shock. Tables 2 and 3 
summarize the results of the spectral analysis of 
MOSl/2 data, using the VNEI model. Spectra 
of some selected regions with the best fit VNEI 
model are displayed in the following figures: Fig- 
ure 5 shows the PN and MOSl/2 spectra of the 
brightest part of the Lobe (region 3), the spectra 
in Figure 6 are extracted from a bright region in 
the shell (region 5), and Figure 7 shows the heav- 
ily absorbed spectra of the CO arm. In all these 
plots, the non-background-subtracted source spec- 
tra are shown, but not the simultaneously modeled 
background spectra. 

The parameters of the VNEI model for the Lobe 
data of the pointing E, with the spectra from the 
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two EPIC detector types PN and MOS fitted sep- 
arately, are consistent with the MOS results of the 
Lobe regions in Table 2 except that the Mg and Si 
abundances from the PN fits turn out to be sub- 
solar in some regions. The most likely explanation 
for this discrepancy between PN and MOS is the 
calibration issue addressed in the EPIC calibration 
status documentation of April 2003 (Kirsch 2003): 
In the range from 0.3 - 1.0 keV, the PN flux is up 
to 10% higher than the MOS flux, whereas above 
1.5 keV, the MOS flux is up to 10% higher. In the 
EPIC spectra, most of the emission lines cannot be 
resolved due to moderate spectral resolution. This 
complicates the determination of the continuum. 
Assuming the same temperature and absorption, 
the flux inconsistency between the MOS and the 
PN can mimic lower abundances for the PN data 
for elements like Mg and Si with lines above ~ 
1 kcV, because the unresolved lines cannot be dis- 
tinguished from the continuum. Although most of 
the fit parameters (absorption, temperature, ion- 
ization timcscale) for the PN, MOSl, and M0S2 
spectra are consistent, the discrepancy between 
the abundances makes us conclude that the re- 
suhs from the MOS data (MOSl and M0S2 data 
fitted simultaneously, if data of both instruments 
are available) are more reliable than the PN re- 
sults in terms of abundances. Therefore, we list 
the MOSl/2 results in Tables 2 and 3. 

3.1. Shell 

3.1.1. Eastern Part of the Shell 

The eastern part of the remnant shell is farthest 
from the pulsar IE 2259-1-586 and has a more cir- 
cular shape. The density of the ambient ISM is 
apparently lower on this side (Kothcs ct al. 2002). 
Thus the eastern portion of the shell seems to be 
the best region to study the least perturbed shock- 
ISM interaction. 

The spectral analysis of regions in the east- 
ern shell (regions 6-9) yields a temperature of 
^ 0.6 keV for the VNEI model, which is slightly 
higher than those in the Lobe and the northern 
portion of the shell, though not significantly. The 
Mg and Si abundances arc consistent with solar, 
and the ionization timescale is ~ 1 x 10^^ s cm~^, 
similar to that in the northern part of the shell 
(see § 3.1.2). The absorbing column density is A^h 
= 5 — 7 X 10^^ cm~^. We also analyze the spec- 



tra of the regions in the eastern shell using the 
model VPSHOCK. The values for A^h , kT, and 
the abundances agree well with the results of the 
VNEI model in all four regions. The lower limit 
of r has been set equal to zero. The upper limit of 
the ionization timescale r is 3.4toi X lO^^ s cm 
for region 6, 4.8^2 5 ^ 10^^ s cm^"^ for region 7, 
3.5^^7° X 10" s cm-3 for region 8, and 6.8t]°i° x 
10^^ s cm~^ for region 9, all higher than and there- 
fore consistent with the average values obtained 
from the VNEI model. 

No indication of nonthermal emission is found 
in these spectra. We model the spectra with an 
additional powcrlaw component and obtain an up- 
per limit for the flux of the powerlaw compo- 
nent. Studies of nonthermal emission from the 
shells of SNRs yield photon indices of around 
2.0: r = 1.8 for Cas A (below 16 keV, Allen 
et al. 1997) and ~ 2.2 for SN 1006 (Bamba et al. 
2003; Allen et al. 2004). Therefore, we assume a 
photon index of F = 2.0 and obtain upper lim- 
its of the nonthermal flux (0.3 10 kcV, 90% 
confidence level): 1.7 x 10^^^ ergs cm^^ s^^ in 
region 6, 1.3 x lO^^'' ergs cm~^ s~^ in region 
7, 3.0 X 10~^^ ergs cm^^ in region 8, and 
1.5 X 10""'^^ ergs cm^^ s^^ in region 9, i.e. <0.2%, 
<1.6%, <0.4%, and <0.2% of the unabsorbed 
VNEI flux, respectively. Gotthelf et al. (2001) 
have studied the emission from the forward shock 
in the northern part of Cas A and report that the 
flux of the nonthermal emission is 50 - 70% of the 
total 0.5 - 10.0 keV flux. In comparison, the upper 
limit of the nonthermal emission from the eastern 
part of the shell of CTB 109 in the energy range 
of 0.5 - 10.0 keV is 0.2 - 1.4%, about two orders 
of magnitude lower than in Cas A. 

In radio continuum, CTB 109 has a very bright 
extended region approximately where the connect- 
ing line from the pulsar through the Lobe inter- 
sects the shell in the northeast (see e.g. Fig. 1 in 
Kothes et al. 2002). In region 5 which lies north 
of this radio bright region, the column density Ah 
(~ 8 X 10^^ cm~^) and the ionization timescale 
T (6 X 10"'^'^ s cm~'^ from the joint fit of PN and 
MOSl/2 data, 4 x 10" s cm'^ from the fit of 
MOSl/2 data, see Table 2) are higher compared to 
the eastern shell regions south of the radio bright 
region, whereas the temperature is significantly 
lower: kT = 0.29 ± 0.03 keV at the 90% level, 
as obtained from both the PN and MOSl/2 fit 
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and the MOSl/2 fit. The kT value for region 
5 differs from that in other regions also at the 
99% confidence level. The spectra are shown in 
Figure 6. Fitting a model similar to other shell 
spectra (iVn = 6 x lO^^ cm'^, kT = 0.6 keV, 
r = 1 X 10^^ s cm~^) results in a modeled spec- 
trum with lower flux below 0.9 kcV and flatter 
continuum above 2 keV, compared to the observed 
spectrum. The optical image of Fesen & Hurford 
(1995) shows that this region covers two bright 
optical filaments with a [Sii]/Ha ratio typical for 
shocked gas rather than of photoionized gas; they 
derive a prc-shock density of 10 20 cm^"*, higher 
than for other optical filaments seen in CTB 109 . 
This might indicate that a shock-cloud interaction 
is occurring in region 5. The ionization timcscalc 
for this region is ~ 4 - 6 times higher than in other 
parts of the remnant shell. Since the regions are all 
located at similar radii, the time that has elapsed 
after these regions have interacted with the shock 
front can be assumed to be comparable. There- 
fore, the higher r in region 5 is consistent with 
the shock having encountered a density enhance- 
ment. 

Although we assume that the undisturbed east- 
ern part of the remnant can be reproduced by a 
model based on the Sedov-Taylor-von Neumann 
similarity solution (Scdov 1959; Taylor 1950; von 
Neumann 1947), the fit results are not satisfac- 
tory if we use the SEDOV model (Borkowski et al. 
2001) in XSPEC. This is presumably due to in- 
homogeneities in the ISM that result in flux and 
spectral variations in the remnant. Therefore, we 
focus on smaller shell segments and analyze arc- 
shaped regions of the eastern shell (see Fig. 3 and 
Table 3) using the evigweight method. As the 
VNEI fit results show, the spectra in the 'outer' 
and 'inner' shell seem to be different: The best fit 
temperature is slightly higher in the outer shell, 
whereas the ionization timescale is lower. The 
foreground absorption is lower for the outer shell 
than for the inner shell. However, if we compare 
the confidence contours of A^h vs. kT and kT vs. 
T for the two shell regions, the 90% confidence 
regions are fully separated, but the 99% regions 
overlap. Thus, there is only a marginal inconsis- 
tency between the temperatures and the ionization 
timescales of the outer shell and the inner shell. 

For a Sedov-phase remnant expanding into a 
homogeneous medium, the temperature and ion- 



ization timescale both increase radially inward. 
This is an idealized case, however. If the kT and 
r differences are real, we may be seeing the effects 
of nonuniformities in the ambient medium. If the 
shock (as seen in projection) encountered a density 
enhancement, that could cause a larger ionization 
timescale and lower temperature, and reduce the 
flux in the inner shell region relative to the outer 
shell region. 

3.1.2. Northern Part of the Shell 

For the northern portion of the shell (regions 
12, 13) Ah « 7 X lO^^ cm'^, kT w 0.5 keV, and 
T PS 1 X 10"'^'^ s cm"''. The temperature seems to 
be lower and the column density A'h higher than 
in the eastern portion of the shell. For regions in- 
side the northern shell segment (regions 11, 14), 
the column density A^h seems to be slightly higher 
(8 — 9 X 10^^ cm~^), perhaps indicating a grad- 
ual increase of the foreground absorption in the 
direction of the CO arm. 

3.1.3. Southern Part of the Shell 

In the southern portion of the shell, there are 
two diffuse spots that appear red in the RGB com- 
posite image (Fig. 2): a bright region in the south- 
east (region 15) and the western tip of the shell in 
the south (region 16). The temperature is 0.6 keV 
and the ionization timescale r ~ 1 x 10^^ s cm~^, 
comparable to the eastern part of the shell. Rather 
than being a result of a lower temperature plasma, 
the red color of these regions appears to be a re- 
sult from a lower Ah (~ 5 x 10^^ cm"^, as de- 
rived from a joint flt of PN and MOSl/2 data, 
~ 5.5 x 10^1 cm-2, from MOSl/2 data only). 

3.2. CO Arm 

We analyze a circular region covering the ex- 
tended tip of the CO arm (region 10). The col- 
umn density Ah is higher than in the eastern part 
of the shefl and the Lobe, A^h > 8 x 10^^ cm~^, 
but the temperature, the Mg and Si abundances, 
and the ionization timescale are comparable to 
the shell of the remnant {kT = 0.5 - 0.6 keV, 
T — 1 ~ 2 X 10^^ s cm~^, as derived from PN and 
MOSl/2 data). 

We also study a non-circular X-ray faint region 
corresponding to the CO arm (see Fig. 3), using 
evigweight. Figure 7 shows the spectra of this 
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region obtained with the EPIC MOSl/2 of the 
pointing PI; the flux falls off more rapidly below 
1 keV compared to other regions of the SNR, con- 
sistent with the presence of additional absorption. 
The result of the spectral fit of the MOSl/2 data 
is given in Table 3; it confirms that the column 
density iVn is almost 10^^ cm~^, higher than in 
other regions of the remnant. The intrinsic spec- 
trum seems to be the same as in most of the other 
regions of the remnant. Therefore, we conclude 
that this region appears fainter because of addi- 
tional absorption. 

3.3. Dark Interior 

There is a low surface brightness region south 
of the Lobe inside the remnant shell. No en- 
hanced CO emission is found in the observations 
by Tatematsu et al. (1990) that coincide with this 
X-ray faint region. We extract the X-ray spec- 
trum of this region using the evigweight method 
excluding the point source which is located within 
this region. The parameters kT, Mg and Si abun- 
dances, and T of the dark region are all consistent 
with the parameters of the CO arm (Table 3), in- 
dicating that the instrinsic spectra of these two 
regions are similar. Only the column density A^h 
is higher in the CO arm. The spectral parame- 
ters also agree well with the parameters of region 
1 (northern tip of the Lobe) which has about the 
same distance to the pulsar as this region. Al- 
though A^H does not differ significantly from that 
in region 1, it seems to be higher in this dark re- 
gion, indicating a higher absorption. 

3.4. Lobe 

There arc three observations which we can use 
to study the Lobe: pointings E, PI, and P2. The 
results from pointings E and P2 are consistent 
with that from pointing PI. The largest number 
of counts is obtained from the observation PI, be- 
cause the exposure time is 2 to 3 times longer than 
in other observations. We also fit all the spectra 
(PN and MOSl/2 data of pointing E, MOSl/2 
data of pointing PI, and M0S2 data of pointing 
P2) simultaneously for each of the regions 1-4. 
The results of these simultaneous fits are all con- 
sistent with the results of the fit of the MOS 1 /2 
data of pointing PI (shown in Table 2). 

In all the Lobe regions, the best fit tempera- 



tures, kT, range from 0.53 to 0.59 keV (from simul- 
taneous fit of PN and MOSl/2 as well as from the 
fit of MOSl/2 data, see Table 2). For Lobe regions 
2 4, the column density TVh « 5 x 10^"'^ cm~^ and 
the ionization timescale t = 2 — 3 x 10^^ s cm~^. 
The rather small differences in temperature, ion- 
ization timescale, and A'h seem to be responsible 
for the apparent spectral differences seen as color 
variations in the RGB composite image (Fig. 2). 
But the detailed spectral analysis indicates that 
there are no statistically significant differences 
among the four regions. In the bright central part 
of the Lobe, Si appears to be slightly overabundant 
(1.6 ± 0.2 times solar), while the Mg abundance is 
consistent with solar. This is not what we would 
expect if the by-products of dust destruction were 
contributing to the X-ray emitting material, be- 
cause the Mg and Si abundances should track each 
other (Itoh 1989; Vancura et al. 1994, and refer- 
ences therein). Since the other three regions have 
Si abundances consistent with solar and all regions 
have Mg abundances consistent with solar, there 
is no compelling evidence for ejecta material in 
these regions. Furthermore, the global NEI fits 
that have all other elemental abundances set to 
solar values produce good fits with no large resid- 
uals around elemental emission lines. 

3.4-1- Line Emission Analysis 

All the EPIC spectra for the Lobe show line 
features which can be attributed to Mg xi, Mg xii, 
and Si xiii lines, whereas the Si XIV line is too faint 
to be visible. Therefore, we use the MgXl triplet 
and Mg xii Lya line as line diagnostics of the NEI 
conditions of the plasma. We analyze the Lobe 
spectra from MOSl/2 data of pointing PI and 
from M0S2 data of pointing P2. The spectrum 
is modeled using a modified 'APEC (Smith et al. 
2001) model that excludes all emission lines^ for 
the continuum and Gaussians for the lines, over 
the energy range of 1.15 - 3.0 keV. The fitted 
parameters for the Gaussians are used to deter- 
mine the line centroids and fluxes. The line en- 
ergies for the Gaussians are modeled as follows: 
1.211 keV for Nex Ly/3, -1.340 keV for Mgxi 
triplet, 1.472 keV for MgXll Lya, ~1.850 keV for 
Si XIII triplet, and ~2.440 keV for Sxv. 

With the spectral resolution of the EPIC de- 
Contact R. Smith for a copy of this model. 
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tectors, wc arc not able to resolve the lines in the 
Mgxi triplet. Therefore, we model the Mgxi line 
feature with one Gaussian and determine the line 
centroid. Variations in the line centroid from re- 
gion to region would be indicative of different rela- 
tive strengths of the resonance, intercombination, 
and forbidden lines. However, the line ccntroids of 
Mg XI are equal for all four Lobe regions within the 
90% confidence intervals (Table 4). The data sim- 
ply arc not of high enough quality to determine 
such a small shift. However, the EPIC spectral 
resolution is sufficient to resolve the Mgxi triplet 
from the Mgxii Lya line. From the fluxes of the 
Mg line complexes, the ratio Mgxi triplet/MgXll 
Lya is calculated. In Table 4, the line energy 
in kcV and the flux in photons/cm^ /s/extraction 
area are given for the Mg xi triplet and the Mg xil 
Lya line, as well as the derived ratio of Mgxi 
triplet/MgXll Lya. 

The resulting lines fluxes and line ratios are 
compared to NEI plasma models obtained using 
neiline (see Appendix A). The calculation pro- 
vides line emissivities for different temperatures 
and ionization timescales. In Figure 8, the MgXl 
triplet/Mgxii Lya line ratio of region 3 for PI 
data derived from the spectral analysis and error 
estimate based on the 90% confidence intervals for 
the line fluxes is plotted as a function of the ioniza- 
tion timescale r and temperature kT. We also plot 
the best fit values from the global VNEI fit for kT 
and T and the 90% errors. Since region 3 has the 
best statistics, it gives the narrowest band in the 
kT-T diagram. Within the errors, the VNEI value 
lies within the allowed region resulting from the 
line ratio analysis. From the diagram, we can de- 
rive the kT range for region 3 as: kT > 0.32 keV. 
This lower limit is consistent with the results from 
the VNEI fit of the EPIC spectra in the total en- 
ergy range of 0.3 - 10.0 keV (Table 2). The ratio 
of the Mg XI triplet to Mg xii Lya docs not provide 
a constraint on the ionization timescale. 

These results show that line diagnostics are pos- 
sible with CCD spectra. Compared to global fits, 
they suffer from fewer systematic uncertainties, 
while it is difficult to quantify the systematic er- 
rors of global fits. Line diagnostics are indepen- 
dent of foreground absorption and offer an addi- 
tional method to constrain the plasma parameters. 
In our analysis the spectral resolution and the pho- 
ton statistics only allow to set a lower limit on 



the temperature. However, the results obtained 
from fitting the lines only and from fitting the en- 
tire EPIC spectrum are consistent with each other. 
This method is a promising technique of analyzing 
plasma conditions and can be applied for brighter 
SNRs as well as for future missions with higher 
spectral resolution and larger collecting area. 

3.4-2. Nonthermal Emission 

In order to check if there is nonthermal emis- 
sion from the Lobe that would be associated with 
a pulsar jet, we fit the spectra of regions 1 - 
4 in the energy range of 0.3 - 10.0 keV with a 
model including a VNEI component and an ad- 
ditional powerlaw component for the Lobe emis- 
sion. The fit does not improve significantly, and 
the flux of the powerlaw component is consistent 
with zero. Therefore, we derive an upper limit 
for the flux of the powerlaw component, assum- 
ing a photon index of F = 1.5; typical photon 
indices for pulsar jets are ~1.2 - 1.7 (Lu et al. 
2002; Gaensler et al. 2002; Pavlov et al. 2003). 
The nonthermal flux of the Lobe (0.3 - 10 kcV) 
is < 9.2 X 10^^^ ergs cm^^ s^^ in region 1 (90% 
confidence level), < 8.0 x 10~^^ ergs cm~^ s~^ in 
region 2, < 7.0 x 10"^''' ergs cm^^ s^^ in region 3, 
and < 6.1 X 10"^^ ergs cm~-^ s~^ in region 4, i.e. 
<0.05%, <0.08%, <0.05%, and <0.07% of the un- 
absorbed VNEI flux in regions 1-4, respectively. 
For the entire Lobe, the upper limit for the non- 
thermal emission is 6.9 x 10~^^ ergs cm~^ s~^, cor- 
responding to 8.2 X 10~^ photons keV""'^ cm~^ s~^ 
at 1 keV. For comparison, the intensity of the Vela 
pulsar jet is ~ 4 x 10~^ photons keV~^ cm~^ s~^ 
at 1 keV (Pavlov et al. 2003). 

4. Estimates Based on the Sedov Solution 

The morphology of CTB 109 and the existence 
of the giant molecular cloud complex in the west 
suggests that the shock wave expanded into a 
dense cloud and decelerated very quickly. To the 
east, the shock wave is expanding into a lower- 
density medium. CTB 109 seems to be confined 
by denser material in the north and the south (Fe- 
scn & Hurford 1995), which might be responsible 
for the smaller radii of the northern and the south- 
ern X-ray shells compared to that of the eastern 
shell (§ 2.1). However, the difference is relatively 
small (~5 - 10%) and overall the X-ray shell of 
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CTB 109 can be considered as a semi-circle. Cal- 
culations for the evolution of SNRs adjacent to a 
molecular cloud show that the shock propagation 
on the low-density side is almost unaffected by the 
presence of the molecular cloud (e.g. Tenorio-Tagle 
et al. 1985; Wang et al. 1992) and can be modeled 
as a spherical problem. The X-ray temperature of 
~3 — SxlO^K and the bright and well-defined 
X-ray rim of the remnant, as well as the lack of 
optical emission except in a few particular regions 
with higher density, indicate that the remnant has 
not yet reached the radiative phase. To estimate 
physical parameters for the remnant, we apply the 
shock jump conditions and the Sedov-Taylor-von 
Neumann similarity solution (Sedov 1959; Taylor 
1950; von Neumann 1947). We assume a distance, 
D, to CTB 109 of 3.0 ± 0.5 kpc as estimated by 
Kothes et al. (2002), and introduce a scaling factor 

~ D/3.0 kpc. For the mass and densities, the 
following relations apply: Assuming cosmic abun- 
dances (Anders & Grevesse 1989), the gas den- 
sity (excluding electrons) is n = 1.1 nn, where 
riH is the atomic H number density. The electron 
number density depends on the ionization state 
of the gas: if we assiimc at least single ioniza- 
tion for each element, Ue = (1.1 — 1.2) uh, and 
the total number density, including electrons, is 
"tot = nc + n = (2.2 — 2.3) nu. The corresponding 
mass per free particle is to = (0.64—0.61) TOp, with 
proton mass TOp = 1.67xl0~^^g. The mean mass 
per nucleus in any ionization stage is Tfin = 1.4 TOp. 
In the following, suffix '0' denotes pre-shock val- 
ues, and suffix 's' denotes post-shock values. 

For the radius of the blast wave, we use the 
distance of the eastern shell from the pulsar as es- 
timated from the EPIC image: Rs = 18'.5 ± l.'O = 
(16 ± 1) fi-i pc = (5.0 ± 0.3) X 10^^ ck cm. If a su- 
pernova occurs near a molecular cloud, the shock 
can be reflected at the cloud surface and cause the 
center of divergence of the velocity field to move 
away from the cloud (Tenorio-Tagle et al. 1985). 
The pulsar might also have a proper motion. Con- 
sequently, in the case of CTB 109 the distance of 
the eastern part of the shell to the pulsar is only an 
estimate of the radius. An alternative way to esti- 
mate the radius of the remnant is to measure the 
half of the extent of the shell in the north-south 
direction. However, as there is dense matter both 
north and south of the remnant that might inter- 
act with the shock front, the north-south extent 



rather provides a lower limit for the size of the 
remnant. 

For the temperature, we use the result from the 
VNEI fits of the spectra in regions 7 to 9. The 
weighted mean of the temperature are calculated 
from the results which have been obtained by fit- 
ting the PN and MOSl/2 spectra separately, and 
from the results of the simultaneous fits of the PN 
and MOSl/2 spectra. As the comparison of these 
results with the values for the eastern shell (Table 
3) show, there is a discrepancy of <0.08 keV be- 
tween the average temperature derived from the 
circular regions and the temperature obtained for 
the outer shell. To allow for the variation between 
the kT estimates, we increase the error estimate 
to 0.10 keV for this calculation, and obtain for the 
temperature of the eastern part of the remnant 
shell: Tx = 0.62 ±0.10 keV = (7.2 ± 1.2) x 10^ K. 

At the discontinuity of the shock front, the den- 
sity jump is n = 4 ng, with uq being the pre- 
shock ambient density of nuclei. Inside the shock 
front, the density decreases towards the center be- 
cause of the adiabatic expansion of the remnant 
(Heiles 1964; Cox & Anderson 1982). In order 
to derive further parameters, we use the distri- 
bution of the normalized density of the nuclei, 
n{R)/no = nB_{R)/nii,o, as calculated by Cox & 
Anderson (1982) and compute the emissivity nu- 
merically. The emitting volume is modeled as a 
cylinder intersecting a spherical remnant with the 
long axis of the cylinder running along the line 
of sight through the remnant. The radius of the 
cylinder is simply the radius of the extraction re- 
gion which, in this case, is r = 80" = 1.2 ^3 pc for 
D = 3.0 ds kpc. 

From the spectral fits of the EPIC data of re- 
gions 7 to 9 with the VNEI model, we derive 
the normalization in XSPEC a.s K = ■j^fjyr x 
/ 1.2 n|_^ dV = (1.6±0.5) x lO'^ cm'^. Since the 
integral of the normalized density n{R)/no over 
the volume can be calculated numerically, the am- 
bient ISM density no is estimated from the nor- 
malization K: 

1.2 x 10-"n?, n f/nuV 
K = ^, "'° / — dV. (1) 

The projected volume is F = 42 d\ pc'^, and 
the integral is /(nn/nH.o)^ dV = 240 dl pc^. 
Thus the pre-shock H density is tih.o = (0.14 ± 
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0.02) (ig cm and the prc-shock density of nu- 
clei no = 1.1 riH.o = (0.16 ±0.02) dg"'^ cm-^. 

In the case of full equilibration between the elec- 
trons and the nuclei, the shock velocity is related 
to the post-shock temperature as: 



T. = 



3m 9 

1! 
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where k = 1.38 x 10"^^ ergs K'^ is the Boltz- 
mann's constant. With a mean mass per free par- 
ticle of m = 0.61 TOp for a fully ionized plasma, the 
shock velocity is estimated as Vg = [(16/cTx)/(3 x 
0.61mp)]5 = 720 ±60 km s-\ 

It is possible that the electrons and ions are not 
fully equilibrated in the shock, in which case, the 
ion temperature Tjon may be higher, and the elec- 
tron temperature much lower than the mean 
plasma temperature T « Tjon- The electrons 
and ions would then equilibrate slowly through 
Coulomb collisions. The plasma X-ray emissiv- 
ity depends primarily on the electron tempera- 
ture, and the deduced velocity would be larger. 
Ghavamian ct al. (2003) and Rakowski et al. 
(2003) have analyzed the blast wave of the SNR 
DEM L71 in the Large Magellanic Cloud and find 
values for Tc/Tjon at different locations ranging 
from 0.01 (almost no equilibration) to 1.0 (full 
equilibration). The evolution of the electron tem- 
perature and its ratio to the mean post-shock 
temperature in two-fluid phase SNR shocks has 
been studied by Itoh (1978). Applying the Se- 
dov similarity solution to the energy equation and 
the equation of state for the electron gas and as- 
suming Coulomb interactions between the elec- 
tron and the ion gas, he shows that the ratio 
g = Te/T is a function of a reduced time vari- 
able u = t3(nQ/£^|;^)^/^^ which describes the ther- 
mal structure of the blast wave; ^3 is the time 
elapsed since the explosion (in units of 10^ yr), 
no the pre-shock density of the nuclei (in units 
of cm^'^), and E^i is the initial blast energy (in 
units of 10^^ ergs). From X-ray measurements one 
obtains the X-ray temperature, Tx, the ambient 
density, no, and the radius, Rs of the SNR. Since 
these three values are related to each other as 
/ = Tx/T = 0.043Tx(i?sno)"^/2^^/^ Itoh (1978) 
derives the value of / = Tx/T from the intersec- 
tion of the curve /{v) with the theoretical curve 
for Coulomb equilibration g{i/) = T^/T. Based on 
this method, we determine a lower limit for the 



ratio Te/Tg « Tx/T^ « 0.4, and an upper limit for 
the velocity = [(16A;Tx)/(3 x 0.4 x O.eimp)]^ = 
1140 ±90 km s-^ 

The age of the remnant can be estimated from 
the shock velocity using the similarity solution: 



(3) 



Under the assumption of full equilibration of elec- 
tron and ion temperatures, this yields an age es- 
timate of t = (2.8 ± 0.3) x 10" ch s = (8.8 ± 
0.9) X lO'' rf-j yr. If equilibration in the shock is 
incomplete, the age estimate would be lower. For 
the lower limit of Tc/Tg ~ 0.4, the age would be 
t = (1.8 ±0.2) X 10" ds s = (5.6 ±0.6) x 10^ dg yr. 

We estimate the initial energy of the explosion 
from 



Rs 



(4) 



where the mean mass of the nuclei is Wn = 1.4 nip. 
Solving for Eq gives Eq = ^^Rsfhnnot~'^ = (7.4± 
2.9) X 10^" ergs, assuming full equilibration 
between the electrons and ions, and Eq = (18.5 ± 
7.3) X 10^° rfg'^ ergs, assuming partial equilibra- 
tion with a lower limit of the temperature ratio of 
Te/Ts = 0.4. The mass swept up by the SNR shock 
wave is M = ^i?^1.4mpno = (97 ± 23) df^ Mq 
if we assume a uniform ambient medium. 

Wang et al. (1992) have numerically modeled 
an SNR next to a dense molecular cloud, using 
the Kompaneets (1960) approximation. They ob- 
tain a geometry comparable to that of CTB 109 , 
1.3 X 10^ yr after the explosion in an interstellar 
medium with a density of no = 0.13 cm~'^. This 
value for the ambient density agrees well with the 
result we obtain from the XMM-Newton data. The 
assumed initial energy is E = 3.6 x 10^" ergs, 2 to 
5 times lower than the energies we derive in our 
calculations. Consequently, the age estimate of 
Wang et al. (1992) is higher. 

Rho & Petre (1997) have fitted a two-temperature 
Raymond & Smith model to the ROSAT PSPC 
spectrum of the northern shell, but the east- 
ern shell is well fitted with a single tempera- 
ture model with a relatively low temperature of 
kT = 0.2lt^;^^ keV. They derive a shock velocity 
of Vs = 380 km s~^, which is much lower than the 
EPIC result, and a remnant age of 1.9 x 10^ yr 
assuming a distance of 4 kpc. This age is 2 - 4 
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times higher than our estimate. If wc assumed a 
larger distance of 4 kpc, the calculation for the 
full equilibration case would yield 1.2 x 10* yr. 

5. Discussion 

5.1. X-ray spectrum 

ROSAT PSPC and Broad Band X-Ray Tele- 
scope (BBXRT) spectra of CTB 109 have been 
analyzed Rho & Petro (1997). Assuming non- 
equilibrium ionization, they derive kT = 1.7 keV 
and n^t = 1.4 x 10^° cm^'^ s for the southern part 
of the remnant. A fit of the Lobe spectrum with 
a two-component Raymond & Smith model yields 
fcTiow = 0.16t°;°!^ keV and fcThigh = OMtoit kcV. 
Their estimate for the column density, Nn, for 
foreground absorption is 9.2tli X 10 cm . No 
spectral lines can be resolved in the ROSATPSPC 
spectrum, whereas the improved spectral resolu- 
tion of XMM-Newton EPIC allows lines to be iden- 
tified, and the continuum can also be better deter- 
mined. Therefore, with the XMM-Newton EPIC 
data we can show that the emission is arising from 
a plasma out of ionization equilibrium and also 
determine the temperature more accurately. In 
contrast to the i^OS'^T results, we obtain a good 
fit of the EPIC spectra of the Lobe with a sin- 
gle NEI model with a temperature comparable to 
that of the high temperature component of Rho 
& Petre (1997). The column density A^h of the 
EPIC fit is lower than that from the i?OS'^T data 
(~ 5 X 10^^ cm~^). Presumably, the lower spectral 
resolution data of i?OS'^TPSPC allowed success- 
ful fits to two equilibrium plasma models, but with 
a higher absorbing column density. 

Parmar et al. (1998) have fitted the BeppoSAX 
LECS spectrum of the whole remnant shell with 
a non-equilibrium ionization model and obtain 
kT = OMto lr keV, Uet = 3.8t? j x 10" cm'^ s, 
and A^H = 6.9t?J x 10^^ cm'^. The analyzed 
spectrum docs not include the emission from the 
pulsar and the Lobe, but encompasses the rest of 
the remnant. The large errors of the fit param- 
eters of the BeppoSAX spectrum are presumably 
caused by the lower spectral resolution compared 
to XMM-Newton EPIC as well as the large extrac- 
tion region including emission from various parts 
of the SNR. 

For the XMM-Newton EPIC data of CTB 109 , 
we use an NEI model for a coUisional plasma with 



variable abundances for Mg and Si, while fixing 
the abundances of other elements to solar values. 
In the shell and in the Lobe, the temperature is 
kT = 0.50 — 0.70 keV, the ionization timcscale 
r « 1 — 3 X 10^^ s cm~^, and the column den- 
sity Ah = 5 - 7 X 10^^ cm~2, with the excep- 
tion of region 5 in the northeastern shell, where 
the temperature is low (~ 0.3 keV) and the ion- 
ization timescale is high (4 — 6 x 10^^ s cm~^). 
In the region corresponding to the CO arm, the 
absorption seems to be higher with A^h = (1-0 ± 
0.1) x 10^^ cm~^, indicating that there is addi- 
tional material along the line of sight to the SNR. 
Abundances are consistent with solar values in the 
whole remnant, except for one region (the bright 
central region of the Lobe), where Si appears to 
be overabundant (~ 1.6 ± 0.2 times solar). Since 
this marginal indication of overabundance is ob- 
tained only for Si in one region within the Lobe, 
while all the other abundances are consistent with 
solar, we conclude that there is no compelling ev- 
idence for cjccta emission. Given the large inten- 
sity variations within the remnant, it is somewhat 
surprising that none of the brighter regions shows 
evidence of enhanced abimdances. It may be the 
case in CTB 109 that the ejecta mass is relatively 
low and the current amount of swept-up matter is 
relatively high such that the cjccta will be difficult 
to detect. Although X-ray emission from ejecta is 
generally observed in SNRs up to an age of a few 
10"^ yr, there arc remnants younger than CTB 109 
which show no evidence for ejecta. For example, 
Kes 79 is estimated to be ~ 6000 yr old and has 
no obvious emission from ejecta (Sun et al. 2004). 

5.2. Interaction Between the SNR and the 
Giant Moleculcir Cloud 

As Tatematsu et al. (1985) have shown, the CO 
emission stretches from the giant molecular cloud 
both north and south of CTB 109. In the opti- 
cal band, filamentary structures have been found 
that are located at the rim of faint Hii regions 
north and southeast of the SNR. At the position 
of the filaments, the H li regions overlap with the 
remnant in projection. This suggests that the op- 
tical emission is mainly caused by the interaction 
of the remnant with dense material harboring the 
Hii regions (Fesen & Hurford 1995). The diffuse 
H II regions are probably associated with the out- 
skirts of the molecular cloud complex to the west. 
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From [Sii] line ratios, Fcscn & Hurford (1995) de- 
rive pre-shock cloud densities of up to ~ 20 cm~^ 
for the optical filaments. New Hi and CO data 
(Kotlics ct al. 2002) make it clear that the rem- 
nant is located in a density gradient: very high 
density in the molecular cloud in the west and low 
density in the emission gap in the east of the SNR. 
As Kothes et al. (2002) point out, the Hi map in- 
dicates that the SNR is not expanding inside a 
stellar wind bubble, since there is no evidence for 
a lower density cavity in Hi with a pronounced 
rim. 

Tenorio-Tagle et al. (1985) have performed two- 
dimensional numerical hydrodynamical calcula- 
tions for the evolution of a supernova remnant in 
or near a molecular cloud. If a supernova with 
an initial energy oi E 10^^ ergs occurs in a low 
density interstellar medium (1 cm~^), the remnant 
shell becomes radiative after 5 x 10^ — 10^ yr 
when it has reached a radius of « 30 pc. How- 
ever, if the explosion occurs in a molecular cloud 
which has a higher density (lO'' cm^''), the Se- 
dov phase is shorter (~ 10"^ yr) and R smaller 
(~ 1 - 6 pc, ShuU 1980; Tenorio-Tagle et al. 1985; 
Chevalier 1999). In the c;ase of a supernova ex- 
plosion inside a molecular cloud, a break-out oc- 
curs if the shock reaches the edge of the cloud. 
For CTB 109 , the position of the pulsar and the 
shell relative to the CMC complex as well as the 
semi-circular shape of the shell indicate that the 
supernova explosion was presumably located out- 
side and close to a giant molecular cloud. As 
Tenorio-Tagle et al. (1985) make clear in their cal- 
culations, a reflected shock forms from the cloud 
surface if the explosion takes place outside the 
cloud, but the shock transmitted into the cloud 
is weak and has a minor effect on the cloud sur- 
face. The reflected shock can cause an additional 
velocity component to the expansion of the SNR, 
directed away from the cloud. Wang et al. (1992) 
have modeled CTB 109 as a supernova explosion 
near (2 pc) a large molecular cloud, assuming an 
initial energy oi E = 3.6 x 10"''° ergs, an inter- 
stellar medium density of no = 0.13 cm~^, and a 
cloud density of n = 36 cm~^. They reproduce a 
semi-circular shell of the observed size at the age 
of 1.3 X lO'' yr. 

The deep EPIC mosaic image of CTB 109 con- 
firms that there is no emission in the western part 
of the remnant. In combination with the morphol- 



ogy of the remnant in the radio, this indicates that 
the shock wave has been stopped completely in the 
west. However, no indication has been found for 
a molecular shock in the CMC (Koralesky et al. 
1998, and references therein). It is possible that 
there is enough material in unshocked parts of the 
CMC which lie in front of the interaction region 
such that any emission from this region could be 
absorbed. Hard point sources discovered in the 
X-ray faint parts of the SNR (§ 2.2) also suggest 
the existence of absorbing material in front of the 
SNR. 

5.3. Lobe as Shocked Cloud 

Our analysis of the XMM-Newton EPIC data 
shows that the spectral properties of the Lobe are 
very similar to those of the remnant shell. More- 
over, the merged XMM-Newton image (Fig. 1) cor- 
roborates that the Lobe and the pulsar are not re- 
lated morphologically, confirming that the Lobe is 
not a jet phenomenon. Rho & Petre (1997) have 
suggested that, due to the interaction between the 
SNR and the CMC in the west, a reflected shock 
has propagated into the SNR. Rayleigh- Taylor in- 
stabilities might have formed between the thin 
SNR plasma and the dense GMC, possibly pro- 
ducing a structure such as the Lobe. Another 
possibility is that there are denser clouds in the 
outer parts of the GMC complex that have been 
shocked by the SNR blast wave. In projection, the 
Lobe and the bright knots in the northern part of 
the remnant seem to trace the outer boundary of 
the CO arm, and therefore might be emission from 
shocked interstellar clouds on the outskirts of the 
GMC complex. 

Tatematsu et al. (1990) point out that a den- 
sity fluctuation in the pre-shock gas might have 
resulted in a local excess of the post-shock emis- 
sion measure seen as the Lobe. They estimate that 
a cloud with a density 2-5 times higher than the 
average pre-shock density in the ambient medium 
could have produced the Lobe. However, they find 
no evidence of CO gas accelerated by the SNR 
shock wave. In the far-infrared, Coe et al. (1989) 
have found a region with bright emission northeast 
of the Lobe. This region also coincides with the 
eastern edge of the CO arm extending from the 
western molecular cloud complex, which suggests 
that this emission is caused by a shocked cloud. 
Fesen & Hurford (1995) have found optical fila- 
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ments at the southwestern edge of the Lobe, i.e. 
closer to the center of the remnant. All the other 
optical filaments detected in their observations are 
located along the northern and southern shell of 
the remnant and seem to be the result of the SNR 
interacting with material associated with the H ii 
regions in the north and the south. Usually, opti- 
cal filaments are believed to be thin, shocked re- 
gions tangent to the line of sight. Since the diffuse 
optical emission in the north extends southward to 
the position of the Lobe, they suggest that the cen- 
tral filamentary emission arises from shocked gas 
along the projected edge of an interstellar cloud. 
Assuming a shock velocity of 100 km s~^, the ob- 
served [Sii] line ratio corresponds to a pre-shock 
cloud density of ridoud = 5 cm~^ (Fesen & Hurford 
1995). 

The interaction between a shock wave and a 
density enhancement like an interstellar cloud has 
been studied theoretically by many authors. Two- 
dimensional hydrodynamical calculations give in- 
formation about the dynamical evolution of an in- 
teraction between the SNR blast wave and a dense 
cloud (e.g. Bcdogni & Woodward 1990; Mac Low 
et al. 1994; Klein et al. 1994). The cloud frag- 
ments due to the interaction. This is also ob- 
served in three-dimensional simulations by Xu & 
Stone (1995), who in addition show that the mor- 
phology of the cloud after the interaction depends 
strongly on the initial shape of the cloud. Klein 
et al. (2003) have performed a laboratory exper- 
iment using the Nova high energy density laser 
at Lawrence Livermore National Laboratory in 
order to study the interaction of a shock wave 
with a high-density sphere located in a low-density 
medium. They confirm that the evolution of the 
sphere after the interaction depends on the density 
ratio of the high-density to low-density medium, 
X, and the Mach number of the shock wave, as 
shown by Klein et al. (1994). The high-density 
sphere (x = 10) is destroyed by the interaction 
with a shock wave with a Mach number of 10. 

In order to better understand the possible ori- 
gin of the Lobe emission in CTB 109 we estimate 
cloud densities from the results of XMM-Newton 
observations based on analytical estimates by Sgro 
(1975) for a planar shock-cloud interaction: There 
are two limiting cases, the 'cold' cloud case and the 
'hot' cloud case, depending on the cooling time 
scale in the post-shock cloud gas relative to the 



crossing time for the transmitted cloud shock. A 
critical density is estimated by equating the cool- 
ing time and the shock crossing time: 

ncrit = 5.8 X 10-4/3^ 4 a- ^ (5) 

where (3 is the ratio of the pressure behind the 
transmitted cloud shock and the pressure in the 
shocked ambient ISM, the shock velocity (in 
km s~^), no the ambient pre-shock density (in 
cm~^), and a the size of the cloud (in pc). For 
a density contrast of x = ^cloud/no = 10, 100, 
and 00, /3 is 2.6, 4.4, and 6.0, respectively (Sgro 
1975). If the cloud is denser than ricrit, the shocked 
cloud gas cools rapidly ('cold' cloud) and, after 
having reached a temperature below ~ 10** K, it 
produces line emission of lower ionization stages of 
the constituent elements. In the 'hot' cloud case, 
the cloud density is lower than ricrit and the tem- 
perature is still higher than ~ 10^ K when the 
transmitted shock reaches the end of the cloud. 
The cloud emits X-rays and optical lines of higher 
ionization stages. However, there will be no region 
that has cooled down to ~ 10^ K, therefore no Ha 
emission will be detected. 

To estimate the properties of a hypothetical 
uniform cloud in the Lobe region, we use the 
values obtained in § 4 from the Sedov analysis 
of the EPIC data, i.e. Vg = 720 km s~^ un- 
der the assumption of full equilibration between 
electrons and ions behind the shock, and no = 
0.16 cm~^. The critical density ncrit as a function 
of the cloud size a (eq. (5)) is plotted in Figure 9 
for three different values for the density contrast 
X = 10, 100, 00. The region below the line corre- 
sponds to the 'hot' cloud case, whereas the region 
above is the 'cold' cloud case. Assuming an am- 
bient density of no = 0.16 cm~^, the initial cloud 
density ndoud = 1-6, 16 for x = 10, 100, respec- 
tively, ncioud is also indicated in the diagram with 
thin horizontal lines. As can be seen in Figure 
9, for X = 100 we would have the 'cold' cloud 
case for all clouds with a > 0.1 pc. However, the 
observed soft X-ray emission from the Lobe indi- 
cates a 'hot' cloud case. If we assume a cloud size 
a = 1 pc and a density contrast x — 10, equation 
(5) yields a critical density of ncrit = 3.7 cm~^, 
while a cloud size of a = 0.5 pc with x = 10 re- 
sults in ncrit = 4.6 cm~^. For a likely cloud size 
a « 1 pc, X = 10 would result in a 'hot' cloud. 
Therefore, the Lobe emission is indicative of an 



16 



interaction between the shock wave of CTB 109 
{vs = 720 km s~^, no = 0.16 cm~'^) and an inter- 
stellar cloud with a density ricioud ^ 5 cm~^. 

6. Summary 

We have studied the Galactic SNR CTB 109 
using EPIC data of five XMM-Newton pointings. 
The deep EPIC image, created from all of these 
observations, shows no emission in the western 
part of the shell. This confirms that the rem- 
nant has its semi-circular morphology because the 
progress of the shock wave has been stopped by 
the CMC. We find no morphological evidence for 
a connection between the Lobe and the pulsar. 
The mosaic RGB image reveals some regions and 
clumps within the remnant which appear to have 
harder or softer spectra than the average spec- 
trum. However, a detailed spectral analysis indi- 
cates that only two regions have significantly dif- 
ferent fitted values for the TVh , temperature, or 
ionization timescale. The Lobe region appears to 
have rather small spectral variations. 

Using the RGB composite image as a guide, we 
extract spectra for different parts of the remnant 
in circular regions as well as polygon regions. Best 
fit results are obtained with a single-temperature 
non-equilibrium ionization model for a coUisional 
plasma. Abundances are determined for Mg and 
Si, while the abundances of the other elements are 
fixed to solar values. As the spectra show no in- 
dication of nonthermal emission, we derive upper 
limits in the energy range of 0.3 - 10.0 keV. In 
the eastern part of the shell where emission is be- 
lieved to arise from the forward shock the upper 
limit of nonthermal emission is less than 2% of the 
thermal emission. 

Analysis of the spectra of four regions in 
the Lobe indicates that the spectral variations 
seen in the RGB composite image are caused 
by small differences in temperature, ionization 
timescale of the plasma, as well as A^h . None 
of these differences are significant at the 90% 
confidence level. There is no significant evi- 
dence for nonthermal emission from the Lobe. 
The upper limits of the nonthermal emission 
are low; the flux upper limit for the entire 
Lobe in the energy range of 0.3 - 10.0 keV is 
6.9 X 10""'^'* ergs cm~^ s"'^ with an average surface 
brightness of 1.1 x 10~^^ ergs cm~^ s~^ arcmin"^. 



about five times lower than the flux from the X- 
ray jet of the Vela pulsar. The detailed spectral 
analysis shows that the spectra of the Lobe and 
the remnant shell arc remarkably similar, with 
kT = 0.50 - 0.70 keV, r = 1 - 3 x 10" s cm'^, 
and A^H = 5 - 7 x 10^^ cm-^. Mg and Si abun- 
dances turn out to be consistent with solar values 
in almost all analyzed regions of the SNR except 
in the bright central part of the Lobe where the Si 
abundance is slightly higher than solar 1.6±0.2 
times solar). This marginal indication of over- 
abundance, only seen in one region of the SNR, is 
not indicative of cjccta emission. The small color 
variations seen in the RGB composite image do 
not result in significant differences in the fitted 
spectral parameters. We use the Mg line emission 
as an NEI line diagnostic. We fit the Mgxi triplet 
and Mgxii Lya line features with Gaussians and 
derive a lower limit for kT which is consistent with 
the results from the global NEI fits of the EPIC 
spectra. 

The detailed spectral analysis docs confirm sig- 
nificant spectral variations for two locations within 
the remnant. The northeastern part of the shell, 
right next to a radio bright knot, contains a bright 
spot with lower temperature {kT « 0.3 keV) 
and higher ionization timescale (r = 4 — 6 x 
10^^ s cm""^). In this region, Fesen & Hurford 
(1995) have found two bright optical filaments. 
The high value for r indicates a region with a pos- 
sibly higher density; it seems that the shock wave 
of the SNR has encountered a particularly dense 
cloud in this part of the shell. The region just 
north of the pulsar corresponding to the CO arm 
has a significantly higher A^h (~ 1.0 x 10^^ cm~^) 
than other regions in the remnant. The tempera- 
ture and ionization timescale arc mostly consistent 
with the rest of the remnant. We conclude that the 
underlying spectrum in this region is not signifi- 
cantly different from the spectrum of the rest of 
the remnant, but the column density of material 
in front of this part of the remnant is higher. We 
note that this larger value of Ah is consistent with 
what Woods et al. (2004) find for their fits to the 
spectrum of IE 2259-1-586 and hence strengthens 
the already-strong case for an association between 
the remnant and the pulsar. In addition, there 
is marginal evidence for a spectral variation with 
radius in the eastern shell of the SNR. The temper- 
ature is higher in the outer eastern shell, whereas 
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the iVu and r values arc higher for the inner sheU. 
This might indicate that the outer part of the shell 
(seen in projection) expands into a medium with 
lower density, while in the inner regions the shock 
is affected by density enhancements in the ISM. 

From the result of the spectral analysis, we 
estimate an SNR blast wave velocity of Vg = 
720 ± 60 km s~^ assuming that the remnant is 
in the Sedov phase and that there is full equili- 
bration of the electron and ion temperatures right 
behind the shock front, and a remnant age of 
t = (8.8 ±0.9) X 10^ da yr, at an assumed distance 
of -D = 3.0 ds kpc. In the case of partial equilibra- 
tion with Tg/Tg = 0.4 as the lower limit, the blast 
wave velocity would be Vg = 1140 ± 90 km s~^ 
and the remnant age t = (5.6 ± 0.6) x 10'^ yr. 
We also calculate a pre-shock density of the nu- 
clei no = (0.16 ± 0.02) cm'^, initial energy 
Eo = (7.4 ±2.9) X 10^° df^ ergs (full equilibration) 
or Eo = (18.5 ±7.3) x 10^° df^ ergs (Te/T^ = 0.4), 
and swept-up mass of M = (97 ± 23) d^-^ Mq. 

The thermal nature of the X-ray emission of 
the Lobe as well as the optical filaments found 
at the southwestern edge of the Lobe are clear 
indications that the Lobe is the result of the 
SNR shock wave encountering an interstellar cloud 
(n-cioud ^ 5 cm~^). We infer from derived shock 
velocities and densities that the shock wave trav- 
eled through the cloud on a timescale which was 
comparable to the cooling time of the shocked gas 
in the cloud. The cloud was heated to temper- 
atures of ^ 10^ K and has not yet cooled down 
substantially, resulting in little optical emission. 

We arc grateful to V. Kaspi for providing 
the XMM-Newton data (Obs. ID 0038140101) of 
IE 2259+586. We would like to thank A. Read 
and M. Freybcrg for useful discussion about the 
XMM-Newton EPIC background analysis. This 
work was supported by the NASA XMM-Newton 
grant NAG5-9914, NASA contracts NAS8-39073 
and NAS8-03060, and in part by NASA grants 
GO0-1127X and GO1-2060X. This research has 
made use of the SIMBAD database, operated at 
CDS, Strasbourg, France. 
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A. Calculating Non-equilibrium Line Emissivities 

neiline calculates line emissivities (using the Raymond & Smith 1977 code in its 1993 update) for certain 
non-equilibrium conditions, specifically an ionizing or recombining plasma with astrophysical abundances. 
We used the code to calculate the line ratio diagnostics described herein, and feel it could be of general 
interest to the SNR community. 

The code is quite simple. After reading in the input values, it initializes the Raymond & Smith (1977: 
1993 update) plasma code, evolves the plasma, and outputs the requested emission. The user can set the 
abundances, the initial electron and ionization temperatures (the initial ionization temperature is used to 
set the initial ionization balance), the pressure, and the maximum time, ionization timescale, or minimum 
temperature for the model. The plasma evolution can be isothermal, isobaric, isochoric, or a special case 
where the plasma is isobaric until the temperature drops to a specified value after which it is isothermal. 
Except in the isothermal case, the electron temperature drops by a user-defined factor (default 2%), and the 
cooling, ionization balance evolution, and line emissivities are calculated by the Raymond & Smith code at 
each step in the evolution. 

The code is written in C and Fortran, and is available at: 
http://cxc.harvard.edU/cont-soft/software/NEIline.l.00.html . 
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Region 3 (east) 
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Fig. 3. — Mosaic of all five pointings towards 
CTB 109 in the energy band of 0.3 - 4.0 keV with 
analyzed regions. For circular regions, numbers as 
used in the paper are given, and the non-circular 
regions analyzed with evigweight are shown with 
dashed lines. A point source is excluded in the re- 
gion 'Dark Interior' (small canceled circle). 




Fig. 5. XMM-Newton EPIC spectra of the 
Lobe region 3 extracted from the pointing E data, 
with VNEI model fits. The position of the Nex, 
Mgxi, and Sixiii lines is indicated. Instead of 
subtracting the background spectrum from the 
source spectrum, it has been modeled simultane- 
ously. The displayed source spectrum and its fit- 
ted model include the background spectrum com- 
ponents. Background spectra are not plotted in 
this and the following spectrum figures. 
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Fig. 1.— See figl.jpg. 



Fig. 2.— See fig2.jpg. 
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Region 1, EPIC MPS 1/2 
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Fig. 4. — XMM-Newton EPIC spectra of the Lobe region 1 extracted from the pointing E data are shown 
with VNEI model fits for PN and MOSl/2. Solid line is used for source + background spectrum, dashed 
line for background spectrum. The background spectrum has not been subtracted from the source spectrum, 
but has been modeled simultaneously and is included in the spectral model of the source spectrum. Strong 
fluorescence lines are visible in the background spectra: Al-K line at 1.5 keV and Cu-Ni-Zn-K line complex 
around 8 keV in the PN spectrum (left), and Al-K at 1.5 keV and Si-K at 1.7 keV in the MOS spectra 
(right). 
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Region 5 (east) 




1.0 2.0 
Energy (keV) 

Fig. 6- XMM-Newton EPIC spectra of the re- 
gion 5 in the eastern part of the sheU extracted 
from the pointing E data, with VNEI model fits. 
For background spectrum see the caption of Fig- 
ure 5. 



Region 5 
0038140101 




10'" 10" to'^ 10'^ 

Tou [s/cm'3] 

Fig. 8. — The hne ratio MgXl /MgXll for region 
3 derived from the spectral fits of PI data is plot- 
ted as function of r and kT as solid line. The 
estimated range for the error is shown with dot- 
ted lines. kT and r values resulting from the 
global spectral analysis (0.3 - 10.0 kcV with VNEI 
model) are plotted with 90% confidence errors. 



CO Arm (0038140101) 
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Fig. 7. - XMM-Newton EPIC MOSl/2 spectra of 
the CO arm (from pointing PI) with VNEI model 
fit. For background spectrum see the caption of 
Figure 5. 




[pc] 



Fig. 9. — The critical density ricrit (Sgro 1975) as 
a function of cloud size a for Vs = 720 km s'~^ 
and no = 0.16 cm~'^ (thick lines). Solid line is 
used for % = 10, dashed line for x = 100, and 
dash-dotted line for x = oo- For x — 10, 100, 
the corresponding cloud density is shown with thin 
lines for no = 0.16 cm~^. The dotted lines indicate 
a = 0.5, 1 and 5 pc. 
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Table 1: XMM-Newton data used for the analysis. 



Obs. ID*" 


Pointing 


direction 


T i b 

Inst." 


Mode'' 


Start time (UT) 


Effective 




RA 


Dec 








Exposure 




(J2000.0) 








[ksec] 


00575401 


23 01 25.2 


58 41 45 


PN 


E. Full 


2002-01-22 17:03:34 


8.0 


(south, S) 






Ml 


Full 




11.0 








M2 


Full 




11.0 


00575402 


23 00 48.8 


59 05 27 


PN 


E. Full 


2002-07-09 07:58:03 


6.0 


(north, N) 






Ml 


Full 




11.0 








M2 


Full 




11.0 


00575403 


23 02 23.1 


58 53 30 


PN 


E. Full 


2002-07-09 13:34:48 


9.0 


(east, E) 






Ml 


Full 




11.0 








M2 


Full 




12.0 


00381401 


23 00 57.6 


58 53 43 


Ml 


Full 


2002-06-11 09:05:01 


34.0 


(PI) 






M2 


Full 




34.0 


01553503 


23 00 57.4 


58 53 37 


M2 


Full 


2002-06-21 09:35:31 


17.0 



(P2) 



"Throughout the paper, the observations are called S, N, E, PI, and P2, as indicated in the table. 

''Instruments; PN; EPIC PN, Ml; EPIC MOSl, M2; EPIC MOS2. All the analyzed data were obtained with the medium filter. 
=Pull: Pull Window Mode, E. Full: Extended Full Window Mode. 



Table 2: XMM-Newton spectral results of circular regions using the spectral model VNEI with solar abun- 
dances except for Mg and Si. Only results from the simultaneous fits of the MOSl and M0S2 data are 
shown, whereas the paper discusses results from PN data as well. Blank sky data are used as background 
data. 



Region 


Pointing 


A^h(10^^) 

[cm~^] 


kT 


Mg 


Si 


nt(lOii) 


xVd.o.f 








[kcV] 


(solar) 


(solar) 


[s cm~3] 






1 


PI 


o 6n+"°2 

"•"'J-o.io 


o.55tEJ:r3 


1 0+0-2 
J^-'J-o.i 


o.9t°:? 


1 4+1-1 


896.1/702 = 


1.3 


2 


PI 


O 

^■^^-0.02 




0.9t- 


1 9+0.2 
-^•^-0.2 


2 0+0-1 
^■"-0.2 


896.6/684 = 


1.3 


3 


PI 


n 40+0.04 


^■"^"3-0.01 


o.9l°:l 


1 6+0-2 


9 r)+0.3 


904.3/721 = 


1.3 


4 


PI 




n fiS+0 03 

"•^•J-o.oi 




1 9+0.2 
^■^-0.2 




896.5/676 = 


1.3 


5 


E 


0.82iro^r 


29+0 03 

^■^''-0.03 


7+0-2 


1 0+0.9 
1-^-0.7 




747.3/600 = 


1.2 


6 


E 


"-'•^■'■-0.09 


0.64lH^ 


1 0+0.4 

^■•^-0.2 


1 0+0-3 




645.9/593 = 


1.1 


7 


E 




0.63l°:l? 


-, 1+0.3 


i.oir4 


0.8l°j 


678.0/552 = 


1.2 


8 


E 


0.67iro^ 


62+0 04 

'-'•0^-0.03 


1 0+0-3 




0.71°:? 


642.8/593 = 


1.1 


9 


E 


0.54t°:^^ 


n 1^7+0-11 

^•^'-0.36 


1 0+0-6 


1 9+0.8 
-•-•^-0.7 




728.7/607 = 


1.2 


10 


PI 


0.88iro^ 


n 'S2+0-01 

"•^^-0.02 




1 0+0-3 
1-^-0.1 


1 7+0.5 
I-' -0.3 


1126.7/918 = 


1.2 


11 


N 


n 70+0.22 


0.531°:°^ 




6+0-2 

U.D_g 2 


^•^-0.3 


637.7/618 = 


1.0 


12 


N 


7n+°-26 
^•'^-0.11 


0.481°:^^ 


1 0+0-1 


1 0+0-3 
-■■•^-0.3 


1 1+2-^ 


653.4/586 = 


1.1 


13 


N 


-0.05 


n ^^2+0-05 




0.9t°| 




733.9/617 = 


1.2 


14 


N 


93+014 










734.4/641 = 


1.1 


15 


S 


0.54lHi 


o.6o;°:°? 


1 0+0-1 


1 0+0-2 




1549.8/1325 = 


= 1.2 


16 


S 


0.57trol 


o.eotro? 


^^•-1-0.1 


n Q+O-3 

u.y_o,2 




1646.6/1475 = 


= 1.1 
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Table 3: XMM-Newton spectral results of arc-shaped region in the eastern part of the shell, the CO arm, 

and the low surface brightness region in the southern interior, using evigweight. The fit results arc from 
the analysis of MOSl/2 data using the VNEl mock^l: abinidancx^s arc fix(xi to solar values except for Mg and 
Si. 'Closed' data arc iisod as background data. 



Region 


Pointing 


[cm~2] 


kT 
[keV] 


Mg 
(solar) 


Si 
(solar) 


nt(lOi^) 
[s cm~^] 






Outer E. Shell 
Inner E. Shell 


E 
E 


u.oa_o 01 

U.0/_o.o6 


0.68iro^ 


-'■•■'--0.1 
■'■•■'--0.2 


o.9tr2 


1-0.2 
1 -1+0.3 
J^-J^-0.3 


566.7/405 = 
347.8/267 = 


1.4 
1.3 


CO Arm 


PI 


0.97tro? 


U.04_0.oi 


u.y_o.2 




1 c+0.4 
-■^■"3-0.2 


918.5/732 = 


1.3 


Dark Interior 


PI 


o.66lS;^? 


u.oo_o.o3 


^■'-'-0.2 




1 c:+0.4 
-■-■3-0.3 


1452.5/1074 = 


= 1.4 



Table 4: Results of the Mg line analysis from XMM-Newton EPIC MOS data (pointings PI and P2) of the 
Lobe. The spectral model is 'APEC No Line' plus Gaussians. 



Region 


Mg 


;xi triplet 


Mj 


5x11 Lya Mgxi triplet/Mgxii Lya 




E [keV] 


Norm*^ 




Norm* Line Ratio 






PI, 


MOSl/2 





1 (North) 


1 QQy+O.OOg 
_o.002 


^■^-0.5 


X 10- 


-4 


2 0+"-'* 


X 10 


-5 


r r+2.5 
"^•3-3.0 


2 (East) 


1 QQO + 0.008 

1.000_Q QQ3 


q 7+1-3 
^- ' -2.0 


X 10" 


-5 


1+°'^ 


X 10- 


-5 


97t~ 


3 (Center) 


1 007+0.004 
1.00 <_0 002 


-■-•^-0.2 


X 10" 


-4 


7 7+5.3 
' •'-4.0 


X 10" 


-6 


161:1' 


4 (Southwest) 


1.0'iU_o 010 


7 9+1.8 
'•"'-2.0 


X 10" 


-5 


i+i-"* 


X 10" 


-5 


79+cx) 
'■^-69 



P2, M0S2 



1 (North) 


1 oqr+0.011 

i.OOO_o.oio 


9Atli 


X 10" 


-5 


1.8+^7 


X 10" 


-5 


c 9+114.8 
■-••^-3.9 


2 (East) 


-1 q^q + 0.013 

l.O1O_0.008 


-r 9+0.2 
-•^•^-0.3 


X 10" 


-4 




X 10- 


-5 


-LO-12 


3 (Center) 


1 qqc+0.006 
i.oo<J_o 005 


-r q+0.5 
-^-•"^-0.2 


X 10" 


~4 




X 10" 


-5 


12 + 168 


4 (Southwest) 


1 q4n+0-0ii 

i.041U_o.oiO 




X 10" 


-5 


o.3t°:«3 


X 10" 


-5 


27tg 



^Norm = total number of photons/cm^ /s/extraction area in the line. 
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This figure "figl.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0408290v2 



This figure "fig2.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/0408290v2 



